A Farewell to Liouvillians 
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We examine the Liouvillian approach to the quantum Hall plateau transition, as introduced 
recently by Sinova, Meden, and Girvin [Phys. Rev. B 62, 2008 (2000)] and developed by Moore, 
Sinova and Zee [Phys. Rev. Lett. 87, 046801 (2001)]. We show that, despite appearances to the 
contrary, the Liouvillian approach is not specific to the quantum mechanics of particles moving in a 
single Landau level: we formulate it for a general disordered single-particle Hamiltonian. We next 
examine the relationship between Liouvillian perturbation theory and conventional calculations of 
disorder-averaged products of Green functions and show that each term in Liouvillian perturbation 
theory corresponds to a specific contribution to the two-particle Green function. As a consequence, 
any Liouvillian approximation scheme may be re-expressed in the language of Green functions. 
We illustrate these ideas by applying Liouvillian methods, including their extension to Nl > 1 
Liouvillian flavors, to random matrix ensembles, using numerical calculations for small integer Nl 
and an analytic analysis for large Nl. We find that behavior at Nl > 1 is different in qualitative 
ways from that at Nl = 1. In particular, the Nl = oo limit expressed using Green functions 
generates a pathological approximation, in which two-particle correlation functions fail to factorize 
correctly at large separations of their energy, and exhibit spurious singularities inside the band of 
random matrix energy levels. We also consider the large Nl treatment of the quantum Hall plateau 
transition, showing that the same undesirable features are present there, too. 



S3 
O 

i o i 

(N 
> 

m 

(N 
(N 

(N 
O 

03 



i 

S3 
O 

o 



X 



I. INTRODUCTION 



It has long been appreciated that the quantum Hall 
effect depends for its existence on localization by dis- 
order of quasiparticles in Landau level tails, and that 
neighboring quantum Hall plateaus are separated by a 
continuous quantum phase transition characterized by a 
diverging localization length.tl |Xhe scaling ideas which 
encapsulate this understandingoa are supported by ex- 
tensive numerical |SUidiesu and by some of the available 
experimental datacm although the full experimental sit- 
uation remains unsettlecB. An analytic theory of the 
transition, however, has proved elusive, even for the sim- 
plest models which include disorder and magnetic field 
but omit electron-electron interactions. 
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In this context, two successive recent developmen 
have attracted interest: Sinova, Meden and G 
(SMG) have introduced a Liouvillian approach to local- 
ization m the lowest Landau level, which Moore, Sinova 
and ZeeB (MSZ) have extended, with the introduction 
of Nl Liouvillian flavors and an, expansion in powers of 
1/Nl about the large Nl limitO. In brief, the Liouvil- 
lian is the time evolution operator for electron probabil- 
ity density, and the information it encodes on localiza- 
tion is integrated over states at all energies within the 
disorder-broadened Landau level. Despite the energy in- 
tegration, the critical behavior of the localization length 
at the plateau transition can, in principle, be extracted 
from the dependence of the Liouvillian propagator on fre- 
quency u>, provided this is known sufficiently accurately 
(see SMG and also Appendix A of this paper) . The orig- 
inal work of SMG reported a calculation of the Liouvil- 
lian propagator using a version of the self-consistent Born 
approximation (SCBA), which is exact for the Nl — oo 



limit of MSZ, and yields diffusive time evolution without 
localization. MSZ have conjectured that the 0(1/Nl) 
correction, which is logarithmically divergent at small u>, 
may be used to arrive at an estimate for the localization 
length exponent v. 

Prima facie there are three reasons to think that this 
work may have outflanked the obstructions faced by pre- 
vious approaches. First, the derivation of the Liouvillian 
formalism by SMG invokes the algebra of density opera- 
tors projected onto the lowest Landau level, thus appear- 
ing to build in the physics of high magnetic fields at the 
first step. Second, the formalism deals directly with a 
disorder-averaged two-particle quantity, avoiding inter- 
mediate calculations of a one-particle Green function, 
and hence, plausibly, goes directly to the heart of the 
matter. Third, the extension to Nl Liouvillian flavors 
is distinct from the standard extension to N orbitals in 
localization problems,£2l which is known to capture weak 
localization pfty-sics but not the quantum Hall derealiza- 
tion transition!!! This gives grounds for hope that the 
logarithm found by MSZ may not simply be due to weak 
localization effects and might indeed be used to get an 
estimate for the correlation length exponent. 

Given the potential importance of these developments, 
it seems useful to investigate the Liouvillian approach in 
some detail and to relate it to established methods. This 
is our objective in the present paper. There are four dis- 
tinct facets to our results. First, we show that although 
the Liouvillian approach has appeared in past work to be 
tailored specifically to the quantum mechanics of parti- 
cles moving in a single Landau level, since it makes use of 
the algebra of projected density operators, it can in fact 
be formulated for any disordered single-particle Hamil- 
tonian. Second, we compare the perturbation expansion 
for the Liouvillian propagator with conventional calcu- 
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lations of the disorder-averaged two-particle Green func- 
tion. We demonstrate that each term in the perturbation 
expansion for the Liouvillian propagator corresponds to 
a specific combination of terms in the Green function ex- 
pansion. As a consequence, any approximation scheme 
within one approach has an equivalent in the other ap- 
proach. Moverover, it is possible to translate between the 
two approaches in either direction: to go from the two- 
particle Green function to the Liouvillian simply involves 
an energy integration, while a more elaborate proceedure, 
which we set out, is required to undo this energy integra- 
tion and pass in the opposite direction. Third, we illus- 
trate these ideas by applying the Liouvillian approach to 
random matrix theory, discussing the Gaussian unitary 
and orthogonal ensembles. We obtain an analytic solu- 
tion in the limit of large flavor number Nl, for arbitrary 
matrix size N, and we supplement this with numerical 
calculations for finite Nl, finding qualitative differences 
between results at Nl = oo (and all Nl > 1) and those 
at Nl — 1. We also undo the energy average at Nl = oo, 
showing that the Liouvillian SCBA has a very different 
character from established approximations when trans- 
lated into a calculation for the two-particle Green func- 
tion. It constitutes an approximation without the usual 
structure based on single-particle self-energies and two- 
particle irreducible vertices. Disappointingly, this is not 
progress as the approximation is pathological: the result- 
ing Green functions exhibit spurious singularities inside 
the band and fail to factorize correctly at large energy 
separations. Finally, we return to the plateau transition. 
Progress in this case is more difficult, because even the 
Liouvillian SCBA of SMG requires a numerical solution. 
We are nevertheless able to show that the two-particle 
Green function generated by the Liouvillian SCBA has 
undesirable features in this case, too, failing to factorise 
as it should for large separations between pairs of its spa- 
tial arguments, and exhibiting spurious singularities as a 
function of energy in this limit. 

In sum, on one hand we have shown generally that any 
Liouvillian approximation can equally well be expressed 
using the better-understood machinery of Green func- 
tions, and on the other hand we have argued that the 
only existing basis for Liouvillian calculations, the X/Nl 
expansion, is seriously flawed. In combination, these re- 
sults leave us pessimistic about the scope for advances in 
the theory of the quantum Hall plateau transition using 
Liouvillian methods. 

In the balance of the paper we set out the technical 
content of these assertions. In Section II we develop Liou- 
villian machinery in a general setting and describe, first, 
how to go from two particle Green functions to Liouvil- 
lians by energy integration, and second, our algorithm 
for undoing this energy integration within a given Liou- 
villian approximation. In Section III we review the 1 /Nl 
expansion scheme introduced by MSZ, again in a general 
setting. Section IV is devoted to a detailed examination 
of the Liouvillian technique applied to the test case of 
random matrix statistics. In Section V we consider the 



quantum Hall problem and describe the pathology that 
is immediately apparent by recourse to the previous re- 
sults. We end with a summary and three appendices 
which provide: a more careful discussion of the critical 
behavior of the quantum Hall Liouvillian than is avail- 
able in previous work; some details omitted in the main 
text; and a construction that yields an algebra for an 
arbitrary single-particle Hamiltonian that is identical to 
the algebra of density operators projected onto the lowest 
Landau level. 



II. GENERAL CONSIDERATIONS 

A. From two-particle Green functions to the 
Liouvillian 

Consider a single-particle Hamiltonian H acting on ba- 
sis states \a) in a space V of dimensionality N . Although 
all of its properties are encapsulated in the corresponding 
one-particle Green functions 

G^E) = (a\ - 4— - 16), (1) 
L — H ±io 

(where S is a positive infinitesimal), in most problems 
of interest their more easily computed disorder averages 
(G^) have little of the interesting information contained 
in G ± . One is then forced to consider higher order corre- 
lators such as the two-particle retarded-advanced Green 
function 

KI;<,b( E >") = G + ab {E + ^)G BA {E - |), (2) 

whose disorder average does contain useful information. 
Analogously one can define K ++ ,K h and K . Here 
and in the following, we use ± as superscripts to indi- 
cate retarded and advanced Green functions, and reserve 
lower case state labels (a, b) for the former, and upper 
case ones (A, B) for the latter. 

The central object in this paper is the energy integral 

/OO j TP 
^KX;bB{E,u). (3) 

We will now show that, quite generally, this can be 
expressed as the one particle Green function of a Liou- 
ville "super-operator" (the Liouvillian) which itself acts 
on the space A\> of all linear operators on the space V. 
Clearly, this Liouvillian Green function is to be distin- 
guished from the usual Green function of the Hamilto- 
nian, introduced in Eq. (|l|). The space Ay is N 2 dimen- 
sional and is spanned by the basis set To empha- 
size that these operators themselves form a linear space, 
we will use the notation \a,b) = |a)(6|. A natural inner 
product on this space is defined as TrP^Q for any oper- 
ators P and Q that belong to Ay. In terms of the basis 
this yields 

(a,b\c,d) = Tr{[\a){b\]^\c){d\} = {a\c){d\b). (4) 
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Finally, pairs of operators P, Q define a super-operator Now consider rewriting Eq. (g) (with tildes denoting 
on A\> via the action PoQ\a, b) = P\a) (b\Q from left and inverse Fourier transforms) 
right respectively. 



/OC j 771 



— r dh I d(2 e'(^/2+«)t 1 + l (B-./2-«s) t2 g+ (ii) g M(t2) 
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i I dte l ^ +lS)t G+ b (t)GB A {-t) = -1 I dte l(w+l8)t {9{t)f{a\e- lHt \b){B\e iHt \A) 



= -i dte t{uJ+l5)t e(t)Tr{(e lHt \A)(a\e- lHt )\b)(B\} = -i dte l ^ +lS)t t (a, A\b, B) (5) 
J -00 Jo 

I 

The single particle Hamiltonian induces a time evolution and the matrix elements of C can be read off from 
super-operator, the Liouvillian C, on Ay 

\a,A) t = e lCt \a,A) (6) 



dt 



\a}{A\=i[H,\a}{A\] = iJ2{H ab 5 BA -5 ab H BA )\b}{B\=iJ2 C aA ., bB \b)(B\ 



(7) 



bB 



bB 



r 



More abstractly, 



C = H ol-loH 



In terms of C we have finally 
Il aA . bB (uj) = (a,A\ 



1 



oj — C + id 



\b,B) 



(8) 



(9) 



which is the promised rewriting of H aA;bB (uj). While 
Eq. (||) does not appear in the work of SMG and MSZ, our 
C is simply the generalisation for arbitrary H of the Liou- 
villian discussed for the quantum Hall plateau transition 
in previous papers. We continue our general discussion 
below and take up the plateau transition in Section |v| 



Liouvillian has at least N zero eigenvalues, since A mm = 
0. Finally, if the eigenvalues of H occupy a band of width 
W, the bandwidth of t is 2W. 



C. Perturbation theory 



From Eq. (M) we can generate the perturbation expan- 



(a,A\C n \b,B) 



(ii) 



B. Exact properties of £ 

A few general and exact properties of C can be inferred 
from the preceding formulae. The eigenvalues X ab of C 
are simply related to the eigenvalues e a of H by 

A Q 6 = e a - e b . (10) 

This implies that X ab = —\ ba , so that the eigenvalues of 
C occur in pairs symmetric about A = 0. In addition, the 



It is interesting to see how this arises from energy inte- 
grating the perturbative expression for the two-particle 
Green function. Consider the contribution to K^. bB {ui) 

at order n in H , which is 

, , n r rrml r rrn — ml 

rn— v 7 x « 7 

(12) 

Integration over E 1 , 
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K 



(n) 

;A:bB 



dE _ 1 ^ 
2~Ki ~~ (oj + iS) n+1 



produces precisely the corresponding term in the expan- 
sion for n(w). 



D. Disorder averaging 

For the Green function perturbation theory with Gaus- 
sian randomness in H , the effect of disorder averaging is 



J 



(-l) m n! 
(n — m)!m! 



[H m ] ab [H n - m ]BA 



(13) 



m=0 
I 



to replace (H 2n ) with a sum of products of all pairwise 
contractions (HH). We can translate this into the lan- 
guage of the Liouvillian theory by using the definition of 
C to express its matrix elements in terms of those of H , 
then using the correlators for the matrix elements of H 
to obtain those for C, 



(£ a A;bB C c C;dJj) — {{^abHsA — H a b5BA)(S c dHDC ~ H cc i5dc)) 

— dabScd(HBAHDc) ~ ^ ab^ DC (H b aH c d) ~ 5BA5cd(H a bHDc) + & 'B A& DC (H abH ~ c d) 

I 



(14) 



and finally using these to average over the powers C 2n 
that occur in the perturbative expansion. 



E. From Liouvillian perturbation theory to Green 
function perturbation theory 

The purpose of this subsection is to construct an al- 
gorithm for passing from a given (likely approximate) 
expression for II(a;) to one for K(E,oj). Specifically, 
we would like to associate uniquely each disorder aver- 
aged diagram for Tl(co) with a corresponding set involving 
Green functions. To make the reader's task simpler we 
we first summarize the basic idea of the algorithm, then 
provide a representative example, and finally state the 
general formula. 



1. Basic idea 

In going between Eq. (|l^) and Eq. ( |l3| ) , we have lost 
the distinction between advanced and retarded propoga- 
tors (since E ± ^ ± iS have been traded for u>) but the 

I 



matrix elements still record where the original disorder 
lines were attached in the two-particle Green function 
diagram. If we reconstruct this information from the Li- 
ouvillian matrix elements, where it must reside since they 
are defined by matrix elements of the Hamiltonian, our 
remaning task is merely one of undoing the signs and 
combinatorial factors introduced h#-.the energy integra- 
tion. This is easier done than said.liil 
2. Example 

Let us see how this works at lowest non-trivial order 
(fourth, since at second order there is only one Liouvillian 
diagram). Of the three diagrams in Fig.|l| we will use 



a B y 

FIG. 1: Diagrams for TI(uj) at fourth order in C. 

the middle one, denoting its contribution to H a A-bB(^) 
by n^. 6B (o;). In terms of H it is (summation over all 
repeated indices is implied) 



n 



(p) 

aA;bB 



(uj) = {CaA;uu£wW;bB){CuU,vv£vV;wW 

[OJ + 10)° 

~ -\5ab{HBwHuA){H 2 



(15) 



IWU 



&BA{HauHwb){H 



(oj + i6) 5 1 

— 2(H au H v t) (HbaH uv ) — 2(H a t,Hvu}(HBvHuA) 

— {H au HBA){H 2 ) u b — (H a bHBu}(H 2 )uA — (H 2 ) au(H u bHBA) — {H a bHuA) (H 2 ) Bl 

+2(H 2 ) ab (H 2 ) B A + 4,(H au HBu)(H ub H UA )] , 



(16) 
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which can be represented graphically by drawing the bare 
Liouvillian propagator as a double line and resolving 
disorder vertices using the Hamiltonian as indicated in 
Fig. |[ This is done in detail in Fig. |[ in which the first 
diagram is Eq. ( |l5| ) and last ten are Eq. ( plf ) . The final 



FIG. 3: Graphical representation of the correspondence be- 
tween Liouvillian and Green function perturbation series. 
First diagram corresponds to Eq. [[5] and contains two Li- 
ouvillian contractions. Next, three diagrams are obtained 
by redrawing one of the contractions using Fig. ||, each of 
which is in turn decomposed by repeating the procedure on 
the remaining Liouvillian contraction. Finally, the last ten 
diagrams represent Eq. FlR 



FIG. 2: Graphical representation oi £. = H o 1 — lo H. 



r *i i. .1 i : i 



step is to undo the signs and combinatorial factors gener- 
ated by the energy integration in going from Eq. ( J12| ) to 
Eq. (|13|). For each term in our example, this may be done 
by simply counting how many times H appears on upper 



FIG. 4: Green functions for the Liouvillian diagram /3 of Fig. 
[l], represented by Eqns. (^6|) and (0). This figure should be 
compared with the last ten diagrams of Fig. pi 



or lower lines, representing retarded or advanced Green 
functions. In this way, writing E± = E±Lu/2±iS, we find 
that Eq. ((]]]) generates a contribution to K^ A . bB (E, u) of 



(/?) (Tr s _ 5 ab {H B wH U A){H 2 )wu , S B A(H au H wz )(H 



■2\ 



aA -' bB( ' ; ~ E\E_ " E+E 5 

— 1 / {H au H v b) (HbaHuv) , (HabHvu)(HBvHjJA) 



4 V E + E - E + E - 

1 ( (HauHBA)(H 2 )ub , (HabHBu)(H 2 )uA , (H 2 ) au(HubH Ba) , {HabHlJ a) {H 2 ) Bl 



4 V E\Et E\E 2 _ E\Et E%E 2 _ 



6 



-2 x - x 
6 



(H 2 



>ab 



'H 2 



' BA . 1 

+4 x - x 

6 



(H au H BU ) (H ub H, 



UAI 



(17) 



This equation is depicted graphically in Fig.||. 

Three aspects of this exercise are worth noting. First, 
a given Liouvillian diagram contains a partial sum of sev- 
eral Green function diagrams. Second, the Green func- 
tion diagrams are all added with positive weights. Third, 
one's conventional intuition about the importance of di- 
agrams in Green function perturbation theory is suspect 
when carried over to Liouvillian perturbation theory. In 
this example, a seemingly simple, non-crossing Liouvil- 
lian diagram actually sums some of the crossing diagrams 
in Green function perturbation theory. 



where Ama are combinatorial coefficients and de- 
notes a particular subset of Wick contractions of (£ 2n ) 
corresponding to \x (and similarly for (} afi in the case 
of (H 2n )). The diagram gives rise to a contribution to 
K 



aA;bB 



{E,uj) of 



y, aM ( 1)m m\(2n-m)l(H m H 



(2n) 



E_ 



+ 



1 j^2n — m+1 



(20) 



The alqorithm 



To apply the method illustrated above to all diagrams 
we need to automate the procedure for keeping track of 
the topology of disorder contractions. This can be accom- 
plished by deforming the problem defined by Eqns. (|^) 
and <M) to the form 



This step can be automated by making the substitutions: 
lu + i8 — > E + E^ 1 p — > xE + /z, and h — > —yE^/z, and 
then, to attach the binomial coefficients, resorting to the 
identity 



C(p,h) = pHol-hloH 
TL{w,p,h) = [u- C(p,h) + iS}- 1 . 



(18) 



Here, p and h are arbitrary parameters which record 
whether H acts on a retarded or an advanced Green func- 
tion line. Consider any particular Liouvillian diagram at 
order 2n, which we label below with /i. Its contribution 
to U(uj) is 



(i 2 



{uj + iS) 



2n+l 



I fjm Tr2n—m\ 

Ea{h) 2n-m h m V± n /"M 
rnaP (u + iS) 2 ^ 1 ' 

mot v ' 



(19) 



dx 



dy 



dz 
2niz 



exp(- 



y + z) 



c m y n 



j. ' n. ! 



(21) 

where the z integration contour encloses the origin once 
in the anticlockwise direction. Thus the algorithm for 
energy un-integrating the Liouvillian perturbation theory 
is reformulated as a particular integral transform. 

To summarize, K(E,u) can be recovered from TI{uj) 
by first generalizing II(ci;) to II(a->,p, h) for a Liouvillian 
as in Eq. (jig) and then carrying out the integrals in 



K{E,w) 



dx 



dy 



dz 
2niz 



exp( 



z)Il(E + E-, -E-t 

z 



-E_) 



(22) 



r 



As the objects on both sides are formally defined by their 
diagrammatic series, this is an exact relationship between 
them. Given a finite set of Liouvillian diagrams, this pro- 
cedure can clearly be carried out diagram by diagram. 
On the other hand, the procedure may prove too cum- 
bersome to deal successfully with a particular approxi- 
mate resummation of Liouvillian perturbation theory. It 
turns out that our program can be carried through for 
the Nl = oo limit of MSZ, applied to random matrix en- 
sembles, as we show in Section [iy|. By contrast, for the 
quantum Hall problem even the Nl = oo limit requires 



numerical solution of an integral equation to obtain the 
Liouvillian. The chore of energy unintegration in this 
case is much more involved and we will consider it only 
in a special limit. It is conceivable that in still other cases 
one may need to resort to numerical resummations of the 
unintegrated series. 
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III. 1/N L EXPANSION 

No exact solution of a problem using the Liouvillian 
formalism is known at present. It is therefore pertinent 
to ask whether there are useful, natural approximations 
in the Liouvillian approach that are hard to uncover in 
the standard approach. The simplest one, suggested by 
the interpretation of the Liouvillian as a random Hamil- 
tonian in its own right, is the self-consistent Born approx- 
imation (SCBA), employed for this purpose by SMG. In 
the supersymmctric functional integral approach of MZS, 
the SCBA is a saddle point approximation, but we will 
instead take the diagrammatic route, in which it is a sum- 
mation of all non-crossing diagrams. In cither case, the 
procedure is justified formally by deforming the prob- 
lem to one with Nl Liouvillian flavors and taking the 
Nl = oo limit. One can then examine the stability of 
the solution by perturbing in 1/Nl, in the hope that 
most of the structure in the problem of interest, Nl = 1, 
survives to large Nl- 

The recipe for introducing Nl flavors into a general Li- 
ouvillan problem is as follows. One replicates the original 
N 2 dimensional space of bilinears Nl times and allows 
for scattering between different flavors. Then 

I& i6fl = (a,A;i\ l \b,B;j), (23) 

u> — L + id 

with 

(£ kl £ mn ) = -^(6 kn S lm +5 km 6 ln )(l kl £ kl ). (24) 
Nl 

where in the Liouvillian correlator we have indicated 
only the flavor indices. The replicated Liouvillian can 
be thought of as an Nl X Nl matrix with individual 
elements which are themselves N 2 x N 2 Liouvillians of 
the original, single flavor, problem constructed from dif- 
ferent realizations of disorder. Moreover, C kl — £ , 
so that (Nl + Nl)/2 independent realizations of disor- 
der are necessary to construct a single realization of the 
replicated problem. The entire NlN 2 x NlN 2 matrix 
is Hermitian. The construction we have described is, in 
the terminology of MSZ, the orthogonal generalization of 
the Nl = 1 problem and is the one they have argued is 
useful in generating a 1/Nl expansion for the localization 
length exponent for the quantum Hall plateau transition. 

To see how the expansion goes, consider as Nl — > oo 
the low order diagrams shown in Fig.^. Evidently, dia- 
grams with crossed disorder lines are suppressed by fac- 
tors of 1/Nl relative to those without crossed lines. This 
suppression of crossed diagrams continues in higher or- 
ders of the expansion. Thus the leading approximation 
at large Nl is to sum all non-crossing terms, which are 
of order N®. As usual, this sum can be carried out by 
solving the equivalent self-consistency equation 

K = SjaSABgi i r ( £ . ft . £) . ft] ij 

aA ' bB U + lS UJ + iS l 1 aA-bB 

(25) 



1 1 1 i 1 




1 i 1 i 1 



FIG. 5: Liouvillian diagrams up to fourth order in C. The last 
diagram contains a crossing and is smaller by a factor 1/Nl 
than the others, which are of order iV£. Only flavor indices 
are shown explicitly. 



FIG. 6: The Liouvillian SCBA in diagrammatic form. Clear 
lines represent bare propagators, and filled ones the full prop- 
agator, II(w). 

represented graphically in Fig.^. At next order (N^ 1 ) 
one needs to sum maximally crossed graphs and rainbow 
diagrams shown in Fig. ||. At any order the Liouvillian 
Green function is diagonal in the flavor index 

^ ;i BM=^H (26) 

and n a A ; bs(w) is taken as an approximation to the Nl = 
1 Liouvillian Green function. 



IV. RANDOM MATRIX THEORY 

We now turn to the problem of random matrix statis- 
tics as a test of the Liouvillian large Nl approach. The 
large number of known exact results on this problem and 
the relative simplicity of calculations make it an ideal 
test case. While there is no phase transition (all Green 
functions are analytic inside the band), we will see that 
there is sufficient non-trivial structure in the statistics of 
the eigenvalues that can be used diagnostically. Oddly 
enough, we will find that the Liouvillian analysis will 
yield some Green functions that are singular inside the 
band (see Section IVE). 

The program is this: for N x N random matrices, we 
introduce an enlarged Liouvillian with Nl flavors. Now 
we have (at least) two different limits to consider. The 
standard one takes N ^S> 1 at Nl = 1, and leads to a 
familiar set of simplifications such the Wigner semicircle 
law for the density of eigenvalues which are summarized 
below. It also leads to a perturbation theory in 1/N 
which is the zero dimensional version of impurity aver- 
aged perturbation theory for single-particle Green func- 
tions in disordered conductors. The new limit we will 
consider is Nl 3> 1 at fixed N. As in all large- n ex- 
pansions, this will prove useful if the large Nl problem 
is sufficiently smoothly connected to the Nl — 1 prob- 
lem of interest. While we will find it useful to then take 
N 3> 1 (but keeping N <C Nl) as well, the possibility of 
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finding non-trivial information on finite matrices is a po- 
tential asset of this limit. Similarly, in the application to 
the quantum Hall plateau transition, the infinite volume 
limit is taken after Nl — * oo. 

Before proceeding we collect some standard definitions 
and results on the random matrix theory for ease of com- 
parison with the following Liouvillian analysis. 



ested in the the density of states (DOS) of H 



p(E) =yS(E- e„) = --ImTr 

t * IT 



i 



E-H + iS_ 



and its disorder averaged correlation function 



(28) 



A. Ensembles 

We will be interested in random Hermitean matrices H 
drawn from the Gaussian unitary ensemble (GUE) if H is 
complex, and the Gaussian orthogonal ensemble (GOE) 
if H is required to be purely real. For these ensembles 
it is sufficient to specifiy the correlators of the matrix 
elements, 



(H ab H cd ) = —(5 a d8bc + aSacSbd) 



(27) 



where a = for GUE and a = 1 for GOE. The variance 
has been normalized by the matrix size, N, to produce a 
bandwidth which remains finite as N — > oo. 



R{E,u) = {p{E + ^)p{E-~)) . (29) 



R(E, to) is related to the two-particle Green functions K 
by 



a A 



(30) 

The Liouvillian Green function allows us to extract the 
Liouvillian density of states (LDOS) 



n(to) = 6(u> — \ r . 



-ImTr 



uj — C + iS 



(31) 



B. Correlators 

In addition to the fundamental one and two particle 
Green functions defined in Section II, we will be inter- 



and it follows from Eq. (|To|) that = / dER(E,uj). 

Disorder averaging (at any N) simplifies the structure 
of various correlators 



{Gmn(E)) — S mn G(E) 
(K aA ;bB(E,uj)) = S ab S A Bl<i(E,uj) + 5 a ASbBK 2 (E,uj) + S aB SAbK 3 (E,uj) 

(n a A;bs(^)) = 8 ab S A B^l(^) + SaASbB^2(^) + 5 aB 5 A bR?,{u) (32) 



with the constraints K2 = -K3, II2 = LT3 for GOE and 
A'3 =113=0 for GUE. Finally, the Liouvillian density 
of states has the form 

= — — IImIIi(w)H — 1(33) 

where the powers of N arise from taking a trace over the 
index structure in Eq. ([32]) . 



C. Standard results: large N at Nl = 1 

At N — 00 the SCBA for the single-particle Green 
function, ■ 

G™ n {E) = 5 mn ^L (l - v/1 - 4«V^ 2 ) (34) 



is exact, and hence the DOS is the Wign er semicircle: 
limj V ^oo(/3( J B))/A r = (l/nv)yjl - {E/2vf. For N large 
but finite, the deviations from this limiting form are 
small. For example, the leading behavior of the DOS 
at energies above the large iV upper band edge, E = 2v, 
is p(E) ~ exp(-4:N(E/v - 2) 3 / 2 /3). Higher order corre- 
lators of the DOS, including R(E,u>), can be discussed 
in each of two limits. First, taking N — > 00 with energy 
arguments fixed and all different, they factorize. For ex- 
ample, with 

Roo{E,lj) = ( Poo (E + |)) ( Poo (E - |)>. (35) 

Second, by scaling the separation of energy arguments 
with the mean level spacing, universal corrleation func- 
tions are obtained, which are dominated at small energy 
separations by level repulsion. Thus, for Nco -C 1, 

R(E,u) oc {uNfp{E) 2 , (36) 
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with = 1 for GOE and (3 = 2 for GUE. 

Combining these with our earlier listing of the gen- 
eral properties of Liouvillians we conclude that the exact 
LDOS will have almost all its weight within the range 
— Av < w < Av for N 3> 1. For = oo it is exactly 
zero outside and vanishes with zero slope (quadratically 
~ (Av ± uj) 2 ) near the edges. Finally, at finite N, there 
should be a pseudogap of width v/N near ui = with 
details depending on the symmetry of the ensemble. 

D. RMT Liouvillian 

1. Nl = oo, N arbitrary 
For the GUE Eq. @ becomes 

10+10 10 + 10 



n2(w) - -]v(^ ni ^ 2 = K^" niH ) 

na(w) = (37) 

The solution is similar to that for N = oo RMT (cf. Eq. 
H). Using Eq. (§|) the LDOS is 

/V 2 _ 1 / ,2 

n( w ) = <sh + — =Ji-_ (38) 

7rV2w 2 V 8v 2 

for |w| < 2v2w and zero otherwise. Notice that the sup- 
port of the spectrum at Nr, — oo is completely indepen- 
dent of N. 

Turning to GOE, the self-consistency condition reduces 
to the quadratic equations 



with solution 



, x 1 2v 2 (N + l), , n9 , , 9 , 

10 + 10 N (LO + 10) 

n 3 H = n 2 n = 2v [ (n 3 H 2 - il» 2 ) = — L- (-1— il» 

iV(cj + idj iv + 1 \c<j + zo 



(39) 



Hi(w) = 
leading to the LDOS 



-Av 2 + (1 + N){N(co + iS) 2 - ^16^ - 8A(1 + N)v 2 {co + iS) 2 + N 2 (co + i<5) 4 ) 

AN(2 + N)v 2 (u> + iS) 



^ 7T V iV+iy n 7 7r(A+l) LO + iS 



(40) 



(41) 



r 



As N — > oo the leading, 0(N 2 ), piece of this result is 
identical to the corresponding piece of 0(u;) for the GUE. 
For finite N, however, the LDOS is zero in the window 
\lo\ < %/2v 2 /A, in stark contrast with the GUE result (or 
anything else one might expect of a disordered system: 
according to this approximation the spectrum of the dis- 
ordered problem is completely gapped) . These results are 
illustrated in Fig. ^. 

There are at least four important differences between 
these exact Nr, = oo results and the corresponding, cor- 
rect, Nr, = 1 solutions: 

• At N L — oo results for both GUE and GOE imply 
a bandwidth of 2-\/2t> for the underlying random 
matrices. This is a factor of v2 less than that of 
the correct solution. 

• The infinite gradient of the Nr, = oo expression 
for Q(oo) at the Liouvillian band edges is at odds 
with the exact result, which has vanishing slope. 



Contributions to the LDOS at the Liouvillian band 
edges come from pairs of random matrix eigenval- 
ues which are near opposite random matrix band 
edges. If one assumes that R will factorize in such 
cases, which is certainly true of the exact result 
and should be expected anyway for such widely 
separated eigenvalues, the form found for LDOS 
implies a DOS divergent near the random matrix 
band edges, with p(E) ~ (u%/2 ± E)' 1 ^. In fact, 
we shall see that R does not factorize at Nr, = oo, 
which turns out to be the central pathology of that 
limit. 

• Although there is a delta function at zero frequency, 
arising from zero eigenvalues of t (see Eq. (|10|)), its 
weight is 1, and not N as it should be. 

• Finally, the two solutions are qualitatively wrong in 
the small frequency limit Qco\ < 1/N) in opposite 
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As is the case with 1/iV corrections to the RMT DOS, 



FIG. 7: LDOS scaled by iV 2 at N L = oo with disorder 
strength v = 1. Solid line is GOE (Eq. [H]) at N = 6 while 
the broken one is GOE (Eq. 0) at N = oo and GUE (Eq. 

at arbitrary N. The delta function at zero frequency has 
been suppressed. 



ways. The GUE LDOS is non-zero with no indica- 
tion of level repulsion, while the GOE develops a 
clean gap over a region of 0(1/N). 

The Nl = oo result above is clearly not an adequate 
approximation to Nl = 1 even for the purposes of com- 
puting properties of the Liouvillian . In the next section 
we will see that matters get worse when we energy unin- 
tegrate the Nl = oo result. Before turning to that, we 
first ask whether one may construct a smooth interpola- 
tion between the two limits by perturbing in 1/Nl- Wc 
offer evidence that this is unlikely, both from considering 
the 1/Nl expansion analytically, and from studying the 
problem numerically for a range of Nl values. 



2. Leading 1/Nl corrections 



As discussed in Section III, the leading correction for 
large Nl to the Nl = oo Liouvillian Green function is 
given by the maximally crossed diagrams. These are re- 
summed by two consecutive geometric series, represented 
graphically in Fig. @. 



FIG. 8: All diagrams for IT(w) at 0{\/N L ). 

Denoting the Nl — oo solution by LIi , the result of 
this calculation for the GUE is 
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FIG. 9: 1/Nl correction to the LDOS for GUE 

the 1/Nl correction to the LDOS is finite near uj = and 
develops a singularity at band edges, as illustrated in Fig. 
|9|. In the case of the GOE, the 1/Nl correction, which 
for brevity we do not display, is divergent at all four 
band edges. These divergences show that any attempt 
to fix the bandwidth problems of the Nl = oo limit, if 
feasible, must involve an analysis including contributions 
at all orders in 1/Nl- We have not pursued this analysis 
further. Instead, we offer evidence from numerical studies 
that expansion about Nl — oo is unlikely to yield useful 
information on behavior at Nl = 1. 



3. Numerics at small Nl 

The simplicity of the random matrix problem affords 
us a different, more direct line of attack via exact diago- 
nalization. Specifically, we have diagonalised numerically 
Liouvillian RMT for a range of Nl and N, both to test 
the relevance of analytic results above and also to search 
for new features, especially in the interesting regime of 
small Nl- Representative plots of Cl(u>) are shown for the 
GUE in Fig.[H] and for the GOE in Fig.||, both cover- 
ing a range in Nl at N = 2. We have searched for, but 
found no significant qualitative differences at larger N. 
Our conclusions from these figures are as follows. First, 
the gross of features of the Nl — oo limit (see Fig. fj]) are 
already observable at Nl as small as Nl = 8. Second, 
Liouvillian theories at Nl > 1 are sufficiently different 
from the theory at Nl — 1, even for Nl — 2 and 3, that 
the quantitative utility of any 1/Nl corrections seems 
questionable. Readers should note especially the rapid 
change in the distribution of small eigenvalues as Nl is 
increased. Third, while one clearly cannot hope to in- 
vestigate analyticity conclusively by these means, there 
are two features to the numerical results which suggest 
behavior non-analytic in 1/Nl at Nl — oo. One of these 
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is the number of states outside the Nl = oo band, which 
appears to decrease exponentially in Nl, as in the band 
tails of the RMT DOS. The other is the presence of oscil- 



lations in the LDOS for the GUE at small energies, as a 
function of Nl: £1(0) vanishes for odd Nl and is non-zero 
for even ones. 



J 



N, =2 



N, =8 




FIG. lOi-i Evolution of fi(w) from Nl = 1 to Nl = 8 for N=2 GUE of random matrices. Solid line is the exact solution for 
Nl = lcJ. Delta functions at zero frequency have been suppressed. 




FIG. lh-i Evolution of Q(uj) from Nl = 1 to Nl = 8 for N=2 GOE of random matrices. Solid line is the exact solution for 
Nl — Delta functions at zero frequency have been suppressed. 



E. Undoing the energy integration 



In this subsection we apply the formalism developed in Sec HE to translate the Nl — oo RMT calculation into 
Green function language. We start from Eq. (|l8|), applied to RMT. Summing the non-crossed diagrams with C(p, h) 
leads for the GUE, to the self-consistency equations, modified from Eq. ( |37|) to 



LL» = 

n 2 (w) = 
n 3 M 



i 



Hx(uj)v 



uj + iS uj + i5 
2phv 2 n 1 (uj) 2 v 2 (p~ hfll 2 {uj) 



N(u + iS) 



0. 



u> + id 



N 



ni(w) 



n 2 (w) 



N 



(42) 



We now take the further limit N — » oo whereon the so- 
lution is 



run 



4(p 2 + h 2 )v 2 



2v 2 (p 2 + h 2 ) 



In this limit the solution for the GOE is identical. 



Using Eq. (|2) we obtain K+^. bB (E, 



5 ab 8ABKi(E,u)) where 
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oo poo 



dx / dyexp(-x-y) 



dze z E^EZ 1 



2-ni2v 2 \(xEZ v ) 2 + (yEZ 1 ) 2 )} 



: 2 - ^{{xEZ 1 ) 2 + (yEZ 1 ) 2 } 



.(43) 



We carry out the integral for large positive E± (see Appendix B for details) and obtain results for general E± by 
analytic continuation. We find 



Ki(E + ,E_) 



1 

2V 



i / E+ J El - Av 2 + £L Je 2 - Av 2 ' 



(44) 



where the signs of the square roots are positive for real E± > 2v, and branch cuts in the complex E+ and E- planes 
run on the real axes, from E± = —2v to E± = 2v. 

The very first thing to note about this result is that, unlike the exact answer 



Ka-mb = G + (E + )G-(E-) = S ab S AB 



E+E- 
Av 4 



Av 2 
E% 



Av 2 



(45) 



it does not factor in its dependence upon E±. At N = oo, 
where the energies E± involve eigenvalues which are at 
infinite separation on the scale of the level spacing, this 
is manifestly wrong. 

To proceed further in this discussion, we need to exam- 
ine the analytic structure of Eq. (f4~4|). Observe that the 
result is an analytic function of real E + and £L when 
both E + > 2 and E- > 2. In this region it is the sum of 
the convergent perturbative expansion in inverse powers 
of E± that is generated by unintegrating the Nl = oo 
series term by term. Also, in this region is purely 

real and the choice of the sign of the imaginary infinites- 
imals in E± is immaterial: it equally well describes an 
approximation to K ++ . 

The series diverges on approaching the boundary of 
this region. In the two-dimensional plane of real E+, 
this happens at the outer boundary of the cross-shape 
formed by the four lines E + — ±2v and E- — ±2i>, 
shown in Fig.|l^. Inside this cross-shape, the choice of 
infinitesimal imaginary parts to E± is essential to spec- 
ify which side one is of the branch cuts in Eq. (f44|). De- 



pending on how this is done, we obtain inside the cross 
both and K ++ , which are no longer equal but are 

both gcnerically complex. The development of imaginary 
parts to the Green functions occurs when one of their ar- 
guments enters the band. Indeed, in the exact expression, 
Eq. (|45|), this is trivially true. 

Random matrices do not exhibit a phase transition 
and so we should expect that K A and K ++ are ana- 
lytic functions of E + and E- except at the band edges, 
E± = Av 2 . This is true of the exact result, Eq. (|45"|). Sur- 
prisingly, it is not true of the correlators derived from 
Eq. (M4). Specifically, inside the cross, while is ana- 
lytic, K ++ exhibits singularities on the circle El + E 2 . = 
Av 2 . For example, along the line E+ = EL, i.e. u> = 0, 
we find 

K+~(E,u = Q) = ^ (46) 

while 



for < E/v < 2 and 



ReK ++ (E, u = 0) = Asgn(- - V2) (47) 
Av v 



^cosbT 1 ( E/v — ] when y/2 < E/v < 2 



ImK++(E,u = 0) = { ' v^CW-v ^ (4g) 

4 sintT 1 ( , E/v ) when < E/v < y/2 



r 



where the negative branch is picked for the cosh . These exact Green functions for comparison, 
results are plotted in Fig. [l3| and Fig. |TJ along with the 
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Finally, one can extract from the knowledge of both 
and K ++ the correlator of the density of states, 
which has the simple form 



R(E,u) 



1 



AltV 2 



e(V - 2E 2 - 



T 



(49) 



It is constant inside the circle E\ 



E 2 _ = iv 2 and zero 
outside. This solves one remaining puzzle. While the Li- 
ouvillian derived two particle Green functions imply the 
same RMT band as the the exact ones, in the sense that 
both are analytic outside the cross, we noted previously 
that the Liouvillian bandwidth was wrong by a factor of 
\/2. We can now recover the LDOS, Eq. ([38]), by inte- 
grating R(E, uj) over energy 



dE 9 uj 2 

Q(2E 2 +—- 

4:TTV 2 2 



1 



~ 8^ 
(50) 



which serves also to verify that our analysis is internally 
consistent. 
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FIG. 12: Analytic structure of K(E^ 
inside the circle. 



E-). R(E,uj) is finite 




FIG. 13: Real part of K ++ (E, E) corresponding to the en- 
ergy resolved Liouvillian solution (solid line) vs. exact result 
(dashed line). 




FIG. 14: Imaginary part of K ++ (E, E) corresponding to 
the energy resolved Liouvillian solution (solid line) vs. exact 
result (dashed line). 



F. Summary of random matrix analysis 

Briefly, in our direct examination of the Liouvillian 
method we find that the large Nl limit is unrepresen- 
tative of the behavior of the Nl — 1 problem. Upon 
energy unintegrating the Liouvillian Nl — oo calcula- 
tion, we find a pathological approximation which fails to 
factor correctly and even produces singularities inside the 
band. Together, these facts indicate that the Liouvillian 
I/Nl expansion is not a useful way of approaching the 
random matrix problem. 

In this analysis we have utilized two complementary 
perspectives on the large Nl method and a comment 
on those is perhaps useful. The first perspective is that 



we generate a family of problems with enlarged symme- 
try groups indexed by an integer each one of which can 
be studied directly — which is what we have done in the 
numerical analysis. The enlarged problems are defined 
by multi-flavor Liouvillians that do not themselves arise 
from from single particle Hamiltonians. The second per- 
spective treats the large Nl limit as formalizing pertur- 
bation theory about a saddle point that is believed to 
capture the relevant physics. In this fashion we generate 
a series in powers of 1 /Nl whose first term comes solely 
from the saddle point. The utility of the method is, of 
course, that setting Nl = 1 in the perturbation theory 
yields the answer for the case of interest. In using the 
Nl = oo answer as input for the energy unintegration 
algorithm we have followed standard practice and simply 
truncated the series at its first term and set Nl = 1. If 
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the large Nl method is useful in describing the Nl = 1 
problem this should be a sensible approximation. Since 
we find that this is not so, the premise must be flawed. 



V. THE QUANTUM HALL TRANSITION 

Our Liouvillian analysis of the random matrix problem 
suggests to us that the method is unreliable and hence 
not to be trusted in its application to the quantum Hall 
problem. This is a disincentive to pursue a similarly de- 
tailed study of the quantum Hall case, especially as it 
would be a substantial undertaking: even the Liouvillian 
Nl = oo analysis of that problem requires a numerical 
solution, and a full energy unintegration would require 
much more intensive numerical work. We therefore con- 
fine ourselves to three sets of observations. First, we 
confirm that our general Liouvillian formalism reduces 
to the one used by SMG and MSZ for the quantum Hall 
case. Second, we exhibit the unintegration algorithm spe- 
cialized to this case. Third, we apply the unintegration 
algorithm to a special sector of the Nl = oo problem, 
that of very large momenta, where it reduces to precisely 
the Liouvillian random matrix problem studied in the 
last section. This then imbeds the pathologies of that 
problem into the quantum Hall case. 



A. Liouvillian formulation 

We consider a charged particle moving in a magnetic 
field on a torus of area 2tt1 2 N, where I = (hc/eB) 1 / 2 
is the magnetic length and N is the integer degeneracy 
of each Landau level. We denote the projection opera- 
tor onto the lowest Landau level by V . As the kinetic 
energy is constant within each Landau level, the Hamil- 
tonian projected to the lowest Landau level is simply 
H = VV{r)V where V(r) is the impurity potential. 

It is convenient to work initially with the states {7 7 |r)} 
which form an overcomplete basis for the lowest Landau 
level; to lighten notation we will write {|r)} for {"P|r)}. 
In terms of these we can define the single-particle Green 
functions in real space 



1 



G ± (E;v 1 ,r 2 ) = (v 1 \ 



1 



(E — H ± iS) 



^1*2} (51) 



and the two-particle retarded-advanced Green function 
in real space 

/f + -(£>;n,r 2 ) = G-{E-u J /2-v 1 ,v 2 )G + (E+uj/2;y^y 1 ) 

where we have chosen a pairing of position coordinates 
appropriate for a diffusion propagator. These are ev- 
idently the analogs of Eqns. (]l|) and (0) in Section II. 
From these we derive the analog of Eq. (g) , 

/OO 
— K + -(E,Lu;r 1 ,v 2 ) 
-oo 2m 



ri,ri 



{uj + iS- C) 



|r 2 ,r 2 ). (53) 



As disorder averaging restores translational invariance, 
we take Fourier transforms and define (denoting the dis- 
order averaged quantities by the same symbols), 



K+-(E,Lj;q 
1 



n(w;q) 



2ttPN 
1 



2irl 2 N 
and the states 



K + ~(E, U) r 1; r 2 )e iq - (r2 - ri) (ir 1 dr 2 
n(w;r 1 ,r 2 )e lq(r2 - ri) rfr 1 dr 2 (54) 



K) = J |r,r)exp(-iq.r)d 2 r , 



(55) 



which represent the projected density operators 
V exp(— i(\x)V in the space of operators A on the lowest 
Landau level. In terms of these we can rewrite Eq. ((52" 
as 



_X+-(£, w;q ) 

-oo 2tti 



— OO 

1 



-<(a 



2ttM 2 N^ Hci1 (lu + iS- C) 



K)>- (56) 



In place of the density operators, it is slightly more 
convenient to use (as a basis for .4) the magnetic trans- 
lation operators r q = e +l 9 / 4 p q , which are closed under 
the algebra 



TqV = exp(zZ 2 q A q'/2)T q+q < 



(57) 



where q A q' — t^qiq^ and are orthogonal with normali- 
sation 



(Tql|Tq 2 ) = iV5 q i, q2 . 

With the identification 



|r q ) 



(58) 



(59) 



one can show that in this basis the matrix elements of 
£ = Hoi-ioHa,re 



i2 qq ' — V^—tj'e 



|q-q'l 2 ! 2 



.qAq'i 2 



= 2iVq-q>e~ 



q-q'l 2 i 2 



sin( 



qAq'l 



e 

>I2 



With these we arrive finally at the form 



n(q,w) 



exp(-g 2 P/2) 
2rf 



1 



lu + iS — C 



(60) 
(61) 



, (62) 



which is the one introduced by SMG (up to a constant 
due to difference in definitons). 
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From known facts about the two particle Green func- 
tion it can be deduced (see Appendix A for details) that 
the exact Liouvillian self energy varies for small q,u> as 



according to diagram's topology. For example, the dia- 
gram of Fig. [l5| is given (aside from energy dependence) 
by 



II -:t (q,u;) = iq 2 D(u)) 



iq 2 LO~< 



(63) 



This is sub-diffusive behavior and is controlled by the 
correlation length exponent as q — > at fixed uj. The 
combined work of SMG and MSZ has yielded a self en- 
ergy, in the 1/Nl expansion, of this form with 



D(uj) ~ D [l 



(64) 



suggestive of an expansion of a critical power law. In the 
following subsection we note that the general algorithm of 
Section II could be used to examine this result at different 
energies in the disorder broadened Landau level. 



B. Energy unintegration 

Formally, the procedure for recovering K (q, E, to) from 
II(q, uj) is a straightforward specialization of Section ||. 
One starts by deforming C to 



d 2 qd 2 p 



(VqV_ q ><VptC p >e- 



1 



-i(qAp+pAQ)r 



and boils down to evaluating a determinant. The wedge 
factors (q A p = z • q x p) appear whenever two momen- 
tum lines cross ( one also selects a sign convention) and 
effectively serve as a fingerprint of each diagram. These 
diagrammatic rules are familiar to workers in the field 
and have enabled considerable progress in series analysis 
of the problemtjlia. Here, they illustrate the connection 




FIG. 15: One of the ten diagrams Fig. |^ redrawn and rela- 
belled using momentum coordinates. 



£qq'(P, h) 



q-q'< 



pe 



he 



This generates a deformation of any given Liouvillian di- 
agram or of some partial summation that leads to an ap- 
proximate form for II (q, u>). In either case, carrying out 
the replacements in Eq. (|l^) and the integrals in Eq. ( p^ ) 
gives the corresponding energy resolved expressions. 

As for the general example of Sec. || the algorithm 
relies on disorder contractions to keep track of the topol- 
ogy of individual diagrams with p and h used to assign 
appropriate weights via Eq. 22. Since the first step of 



the procedure is particularly transparent in the LL set- 
ting we briefly consider the analog of Fig. ||. In par- 
ticular, for each of the ten "unfolded" diagrams we can 
immdcdiatcly write down its value by drawing momen- 
tum carrying lines (both disorder and external) inside the 
"particle-hole circle" and integrating correlators of disor- 
der potential with an exponential/phase factor chosen 



between Liouvillian and Green function perturbation se- 
ries. 

In order to apply this procedure in full to the Liouvil- 
lian Nl = oo limit, we must compute II(q, oj,p, h), which 
in turn requires numerical solution of an integral equation 
and is therefore a substantial undertaking. We turn now 
to the large q limit, for which the energy unintegration 
can be carried out more easily. 



C. Energy unintegration at Nl = oo, as q — > oo 

The key observation! that allows us to execute the 
energy unintegration procedure in this limit is that the 
self-consistency equation simplifies drastically. Defining 
II(q, uj) by II(q, u) = e~ l q / 2 II(q, u), as q — > oo the self- 
consistency equation 



J 



n(q,cj) = 



i | n(q,^) 

uj + iS uj + iS 



4^2 N^q.q+q ■'-q+q' 



, q )II(q+q',w) 



(66) 



reduces to 



n(oo, uu) 



1 



2w 2 n (oo,cj) 



uj + iS lj + iS ' 

which is identical to the corresponding RMT equation (ei- 
ther GUE or GOE) at N = oo. It follows then that the 



r 



(67) g lven 



energy unintegrated solution of this equation, the quan- 
tity i?(co, £+,£_) = lim q ^ 0O e+9 2 ' 2 / 2 X(q,^+,S_), is 



i (OO, = limq^oo e H A (C[, ti+ , 1 

by the expression already exhibited in Eq. ( 
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K(oo,E+,E_) = -L 



sm 



El 



4v 2 +E-JE 2 _ -4v 2 



E 2 + 



E 2 



4v 2 



TV 

2 



(68) 



As before, the two particle Green function fails to fac- also correspond to large spatial separations in magnetic 
torize in energy. This has, however, a more serious con- fields. More precisely, one has 
sequence in the quantum Hall problem. Large momenta 



and hence 



(G + (t,t,E+)G-(<0,0,E-)) = J ^exp(-^\zxr + ql\ 2 )K(q,E + ,E_), (69) 



lim (G(r, r, E+)G(0, 0, £_)) = lim ^K(q, E + ,E_) . (70) 

r^oo q-+oo /Til 



So we discover that lack of energy factorization translates 
into a lack of factorization at infinite spatial separation, 
which is clearly unphysical. In addition, we again find 
that the retarded-retarded Green function is singular in- 
side the impurity band in a limit where the exact answer 
is neccessarily analytic. 



VI. SUMMARY AND OUTLOOK 

As this has been a largely technical discussion of the 
Liouvillian formalism, it is perhaps useful to summarize 
the main argument again. The work of SMG and MSZ 
has suggested that the critical divergence of the localiza- 
tion length for the (non-interacting) quantum Hall tran- 
sition can be computed within the Liouvillian formalism 
by the 1/Nl expansion in the number of Liouvillian fla- 
vors. We find, based on testing the Liouvillian approach 
on random martrix theory and energy unintegrating it at 
large momenta in the quantum Hall problem, that this 
program has serious problems already at Nl = oo which 
do not appear to be remediable within the 1 /Nl expan- 
sion. Accordingly, we conclude that the calculation of 
MSZ does not represent a computation of the quantum 
Hall correlation length exponent. As there isn't another 
approximation scheme that readily suggests itself in the 
Liouvillian formalism, we are pessimistic, as indicated 
by our title, about its utility for making progress on the 
problem of the quantum Hall transition. 

That said, there are a couple of loose ends in our analy- 
sis that would be nice to tie up. First, we have not carried 
out the energy unintegration of the Liouvillian Nl = oo 
answer for the quantum Hall problem near q = and it 
would be of some technical interest to see if that has an 



analytic structure different from large momenta as well 
as a technical challenge to see how that might get done. 
Second, it may be possible to use our ideas on energy 
unintegration to see if the logarithm found by MSZ has 
any interpretation in Green function perturbation the- 
ory. Perhaps an intrepid reader will be inspired to sort 
these out. 

Finally,_as we were finishing this work, there appeared 
a preprinto by Moore which studies the evolution of the 
QH Liouvillian as a function of Nl. Its most relevant 
finding, via numerical analysis at small Nl, is that the 
Liouvillian theory likely exhibits metallic diffusion, in- 
stead of critical scaling, at all Nl > 1. Evidentl y, this 
is consistent with our own conclusion in Section IV D 3 



that the theory for all Nl^I is unrepresentative of the 
original problem of interest (Nl = 1). 
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APPENDIX A: CRITICAL BEHAVIOR OF THE 
QUANTUM HALL LIOUVILLIAN 

For completeness, but also because some of the original 
discussionQ was incorrecto, we discuss here how it is that 
the energy integrated Liouvillian propagator can be used 
to extract the localization length exponent. One might 
worry that, as the critical wavefunctions are a set of mea- 
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sure zero in the spectrum, any information on them will 
be lost upon integrating over the band and the behavior 
at small ui and q will be that of an insulator. This is not 
the case: the Liouvillian Green function exhibits critical 
asymptotics intermediate between those of the metal and 
insulator in the transport limit q <C to, with the diffusion 
constant D(u>, q_^ 0) ~ cj 1 / 2 ", where v is the correlation 
length exponento. 

Near the critical energy E c , the two particle 
Green function Jbr the quantum Hall problem can be 
parametrized aala 



K(E,u;q) = 



2irp(E) 



-D(E,Lo;q)q 2 



(Al) 



In the transport limit Dq 2 -C u> we can expand the de- 
nominator to obtain 



K(E,u;q)^2irp(E) 



1 Dq 2 



These forms also hold for a metal, with D(E,co;q = 0) 
non-zero, and for an insulator with D(E,u>;q = 0) van- 
ishing for u> — > roughly as u> 2 . For the quantum Hall 
problem, D(E,u,q) has a non-trivial dependence on E, 
lo and q. 



To quantify this dependence, consider the two impor- 
tant length scales in the transport limit. These are 
the localization length £ ~ l\E c /(E — E c )\ v and L u = 
{p{e)oj)~ 1 ^ 2 which is the size of a box with level spacing 
of order the frequency. Then, when L u <C £(E), localiza- 
tion is unimportant and we may take D{E, to, q) rs D Cl a 
non-zero constant. In the opposite limit, L u 3> £(E), we 
have the insulating behavior, D(E, u>, q) = D ins (uj) ~ ui 2 . 
In the following we will assume a sharp crossover between 
these two limiting forms, which is sufficient for our pur- 
poses. 

Armed with these facts we can now carry out the en- 
ergy integration to deduce the form of n(q, ui) in the 
transport limit. Specifically, 



(A2) J ' dEBeK(E,u;q) = 2irl s J dEp(E)D{E, u, q = 0) 



(A3) 

We will see that the integral is dominated by energies 
near E c so that we can replace p(E) w p(E c ) and trade 
the integral over E for one over the localization length 
using the asymptotic relation £ ~ l\E c /(E — E c )\ v . This 
yields 



'l rp ll/v poo 

ImH(q,c) ~ p(E c ) / d£C 1 - X/v D{E c + E c (l/0 1/u , «, q = 0) 

LU 2 V J l 



1 



EJ 1 / 1 



,-,2 



P(E C ) 



i Jl^ 



^E c p(E c ) Ans(w) + D c {l'p{E c )u) 



l/2v 



D c q 2 



or 



(PpiE^u;) 1 ^ 



r 



(A4) 



where the last simplifcation follows upon noting that 
Z?i ns (w) ~ to 2 is much smaller than u>~ for v > 1/4, 
which is true in general of a random critical point in two 
dimensions on the grounds of the Harris criterion, and 
holds for the reasonably precise estimates of v available 



J 



in this case. As a statement about Liouvillian theory, 
this result implies that the self-energy of the exact Liou- 
viliian Green function has the behavior q — > 0) = 
ui — 1/H(u),q) ~ iq 2 LO^ . 



APPENDIX B: DETAILS OF ENERGY UNINTEGRATION 



We start from Eq. ([43|). After a change of variables, it reads 

rdr exp(— r(£J_ cos 9 + E+ sin 9)) f dze z 



K{E +1 E_) = I d9 

Jo Jo 



d9 

ii JO 



1 

2^ 



2v 2 r 2 

dr exp(— r(£L cos 9 + E + sin 9)) 

v 

Rcos(9 + d>) \ 



2ni 
h{2vr) 



[z - ^z 2 - (2vr) 2 } 



d9 

\^/R 2 cos 2 ( 



(Bl) 
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where we have defined 



and 



R = ^\] E ++ E - = ^ VE 2 + (cj/2 + io+y 



tan i 



El 



(B2) 



(B3) 



Proceeding further, 



K(E + ,E_) 



which is Eq. (| 



1 

2^ 

1 

1 

2^2 



_ 1 R cos c 



sin 



j Rsmcj) 7r 
VR 2 - 1 ~ 2 



sin 



E 2 - 4v 2 



E+\ Ei. — 4« 2 + E-\/E_ -Av 



Ei+Ei- 4v 2 



E 2 - 4v 2 



(B4) 



APPENDIX C: MAGNETIC ALGEBRA AND 
RANDOM MATRICES 



so that 



In this appendix we demonstrate that an algebra like 
that of density operators in the lowest Landau level may 
be set up for an arbitrary Hamiltonian. We consider 
an orthonormal basis {\a}} for V and the corresponding 
basis {\a, b)} on A\>, where a and b range from 1 to N). 
We define 

|q) = j^j2 eW/Ar E + *v> ny q * n2 * /N (ci) 



for integer q x and q y (ranging from 1 — > N) which also 
form an orthonormal basis. We can rewrite 



H = ^2H mn \m,n) = ^V_ q |q), 



(C2) 



idt 



N l/2 

P 



v- 2 , y _ sin ^AP| p + q) . 



N 



(C5) 



These expressions parallel those for the Liouvillian in the 
lowest Landau level. Strikingly we have recast a generic 
problem in the notation of the magnetic translation oper- 
ators. Of course, the definition of the (discrete) momenta 
is adhoc, and so, while this rewriting appears suggestive, 
it is not generally useful. Nevertheless, the GUE retains 
its simplicity in this approach: Eq. (27) implies that 



where 



V- 



H, 



/VI 



+<?«■« -iqx{n+q v /2)2ir/N 



(C3) 



<W = -<5 p ,_ q . 



(C6) 



Switching to operator notation, with r q = |q), one can 
check that 



[Tq,r p ] = 2isin 



7rq A p 

N 



'q+p 



(C4) 



Disorder-correlators are thus diagonal in q, continuing 
the analogy with spatially extended systems. 
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Note that Eq. ( p7j ) is equivalent, for general N, 
to matrices being distributed according to V[H] ~ 



exp[- 



N 



(1+a) 



^Tr{H 2 }}. 
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